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We investigate magnetically-induced ferroelectricity in Cu2MnSnS4 by means of Landau theory of 
phase transitions and of ah initio density functional theory. As expected from the Landau approach, 
ah initio calculations show that a non-zero ferroelectric polarization P along the y direction (of the 
order of a tenth of /xC/cm 2 ) is induced by the peculiar antiferromagnetic configuration of Mn spins 
occurring in Cu2MnSnS4- The comparison between P, calculated either via density- functional- 
theory or according to Landau approach, clearly shows that ferroelectricity is mainly driven by 
Heisenberg-exchange terms and only to a minor extent by relativist ic terms. At variance with 
previous examples of collinear antiferromagnets with magnetically-induced ferroelectricity (such as 
AFM-E HoMnOs), the ionic displacements occurring upon magnetic ordering are very small, so that 
the exchange- striction mechanism (i.e. displacement of ions so as to minimize the magnetic coupling 
energy) is not effective here. Rather, the microscopic mechanism at the basis of polarization has 
mostly an electronic origin. In this framework, we propose the small magnetic moment at Cu sites 
induced by neighboring Mn magnetic moments to play a relevant role in inducing P. Finally, we 
investigate the effect of the anion by comparing Cu2MnSnSe4 and Cu2MnSnS4: Se-4p states, more 
delocalized compared to S-3p states, are able to better mediate the Mn-Mn interaction, in turn 
leading to a higher ferroelectric polarization in the Se-based compound. 

PACS numbers: 75.50.E, 77.80.-e, 75.85.+t 



I. INTRODUCTION 

Multiferroic materials have attracted both scientific 
and industrial interests, due to their profound physics 
and novel functionalities, respectively. |lj Among them, 
many studies have been recently performed for "im- 
proper" ferroelectric materials, [2|46| where ferroelectric- 
ity is driven by non-centrosymmetric spin, charge, or 
orbital ordering. In this case, a stronger magnetoelec- 
tric coupling is expected, compared to the conventional 
proper ferroelectrics. So far, mainly manganites (such as 
TbMn0 3 , TbMn 2 5 , and HoMnOsB-Q) and other per- 
ovskite oxides 0, H[ have been carefully investigated by 
means of theoretical and experimental approaches, in or- 
der to better understand cross-coupling effects between 
spin/orbital orderings and dielectric properties. 

The search for novel multiferroics other than oxides has 
already started: recently Nenert et a/.[9[ have suggested 
the ternary copper chalcogenide Cu2MnSnS4 (with nom- 
inal valences as Cu 1+ , Mn 2+ , Sn 4+ , and S 2_ ) as a can- 
didate of "improper" ferroelectricity, using a symmetry- 
based analysis and the Landau phenomenological theory, 
with the ferroelectricity induced by a peculiar antiferro- 
magnetic (AFM) Mn spin configuration. Experimentally, 
neutron diffraction measurements have clarified that such 
AFM configuration occurs under T/v=8.8 K, with Mn- 
spins slightly deviating from the c axis by an angle be- 
tween 6° and 16° degree. Moreover, the magnetization 
curve has suggested the spin quantization axis to be close 
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to the c direction. |10[ Additionally, the small dependence 
of susceptibility upon the direction of magnetic field im- 
plied a rather small magnetic anisotropy. 

Similar to the famous stannite, Cu9FeSnSzi,[ll| 
Cu2MnSnS4 crystallizes in the tetragonal sphalerite su- 
perstructure, which shows the /42m (No. 121) space 
group with symmetry operations {E, C^x, ^2^, C^z, 
&dx, &dy, S4, ^ 1 }.[l2| The crystal structure is similar 
to zincblende, often seen in III-V and II- VI type semi- 
conductors, where each cation is tetrahedrally coordi- 
nated to anions. In this respect, multiferroic effects in 
this class of materials might prove useful in the context 
of semiconductor-based spintronics. Although the spha- 
lerite crystal structure itself lacks inversion symmetry, 
this does not automatically imply that the crystal is po- 
lar: due to the mirror and 7r-rotation symmetries, the 
ferroelectric polarization P is forbidden along any direc- 
tion. On the other hand, Nenert et at. have claimed that 
the experimentally observed AFM spin order breaks some 
symmetries, so that P arises along the y direction (i.e. 
the crystallographic b direction) with possibly large mag- 
netoelectric coupling. [9] The Landau phenomenological 
theory, based on the symmetry analysis of both the crys- 
tal structure and the magnetic order, is a powerful tool 
to investigate magnetically-induced ferroelectricity. (El fl3| 
However, in order to quantify the value of P and to clar- 
ify the relation between magnetism and ferroelectricity 
via the electronic structure, a combination of Landau 
theory and ab initio density functional theory (DFT) 
calculations is desirable. [4[ In this paper, we investigate 
magnetically-induced ferroelectricity in Cu2MnSnS4 by 
using both Landau theory and DFT in the aim of care- 
fully understanding the microscopic mechanism at the 
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FIG. 1: (a) Nuclear crystal structure of Cu2MnSnS4- (b) and 
(c) indicate AFM unit cells with the propagation vector k — 
(§, 0, |) corresponding to ferroelectric domains with positive 
P and negative -P, respectively. For simplicity, only Mn 
atoms are shown in (b) and (c). Orange arrows at Mn sites 
indicate the spin directions. The magnetic unit cell is doubled 
along the a and c directions, compared to the nuclear unit cell. 
Labels of Mn spins are shown in each sphere. 



II. METHODOLOGY AND STRUCTURAL 
DETAILS 

DFT calculations are performed for Cu2MnSnS4 and 
Cu2MnSnSe4 by using the "Vienna Ab initio Simula- 
tion Package" (VASP), based on Projector Augmented 
Wave (PAW) pseudopotentials.Q The Perdew-Becke- 
Erzenhof (PBE) approach to the generalized gradient 
approximation (GGA) is employed for the exchange- 
correlation potential. [15| The plane wave energy cut-off is 
400 eV for each atom. During structural optimization, a 
threshold on the atomic forces is set as 0.01 eV/A. Inter- 
nal atomic coordinates are optimized starting from exper- 
imental data of Cu 2 MnSnS 4 [Io| and Cu 2 MnSnSe 4 .[l6| 
The 2x4x1 Monkhost-Pack /c-point grid in the Brillouin 
zone is used.[l7| Cu-3<i4s, Mn-3<i4s, Sn-5s5p, S-3s3p, 
and Se-4s4p electrons are treated as valence states. The 
GGA+£7 calculations within the Dudarev's approach fl8| 
are performed by applying a Hubbard-like potential only 
for Mn-d states. An on-site Coulomb parameter U = 4.0 
eV and an exchange parameter J = 0.89 eV are used. 

Lattice constants are fixed to the experimental data: 
a = b = 5.514 A, c = 10.789 A for Cu 2 MnSnS 4 [Io| and 
a = b = 5.736 A, c_ = 11.401 A for Cu 2 MnSnSe 4 [3 
in the tetragonal 742m structure. In order to im- 
pose the experimentally observed AFM spin configu- 
ration in Cu 2 MnSnS 4 [10] with the propagation vector 
k = (^,0, -|), we build the magnetic super-cell with 8 
formula/units (i.e. 64 atoms/cell, see FigOJb) an d (c)) 
(note that, to our knowledge, the experimental magnetic 
configuration for Cu 2 MnSnSe 4 has not been reported 
yet). Since the observed AFM configuration breaks the 
symmetries (except for C 2y ), ferroelectric polarization is 
allowed to be magnetically-induced along the y direction. 
To evaluate P, we also perform calculations for the ferro- 
magnetic (FM) spin configuration, which is considered as 
a reference non-polar state: the ferroelectric polarization 
is then calculated as AP = Pafm — Pfm- To connect 
paraelectric (PE) and ferroelectric (FE) states, calcula- 
tions for noncollinear spin configurations (see below) are 
also performed according to Ref.pjl]. The Berry phase 
approach developed by King-Smith and Vanderbilt |20| is 
employed to calculate the electric polarization P, where 
we integrate over eight /c-point strings parallel to P. 



basis of the multiferroic behaviour. To check the effect 
of the anions on P, we perform similar simulations by 
substituting S with Se, i.e. for Cu 2 MnSnSe 4 . First, we 
will discuss the direction and the value of P obtained by 
DFT, confirming the predictions based on Landau the- 
ory; afterwards, we will focus on the ferroelectric switch- 
ing, in terms of different Mn-spin configurations along the 
adiabatic path between positive and negative ferroelec- 
tric states. Finally, the microscopic mechanism leading 
to ferroelectricity via spatially-polar Cu and Mn charge 
densities caused by Heisenberg-exchange is discussed. 



III. MAGNETICALLY INDUCED 
FERROELECTRIC POLARIZATION 

A. Landau theory 

In the Landau potential terms, the ferroelectric polar- 
ization P appears as coupled with magnetic order pa- 
rameters, in terms which are invariant under symmetry 
operations. [2lL l22| Nenert et aZ.[j| have built up the free 
energy equation by considering two Mn spins (S\ with 
position (0,0,0) and £ 2 with position (\,\,\) in the nu- 
clear unit cell) and set up two magnetic order param- 
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eters, e.g., the FM order parameter, M = Si + S2, 
and the AFM order parameter, L = Si — £2- From 
the free energy, they have showed that only the P y com- 
ponent is coupled with both L x and L z , and also with 
My. Thereby, they have concluded that a spontaneous 
polarization arises along the y direction upon the AFM 
configuration. Here we extend their model to fit into our 
AFM super cell with eight Mn spins, so that we can dis- 
cuss the adiabatic switching path between the positive 
FE (+.P) and negative FE (— P) states by changing the 
direction of Mn spins and the consequent electric polar- 
ization. (The same procedure was applied for HoMnOs 
with E-type AFM configuration. 0, [H|) With these eight 
Mn spins (labelled as Si, • • • , S 8 , as in FigH]), two AFM 
order parameters are constructed as 

Ei = Si — S2 — S3 + £4 — £5 + #6 — £7 + S$, 
E2 = Si — S2 — S3 + £4 + £5 — Sq + £7 — S$, (1) 

corresponding to the positive and negative ferroelectric 
domain phases with +P and — P, respectively. In order 
to switch between the two phases (as seen in the fol- 
lowing section), four body-centered Mn spins (£5, Sq, 
£7, Sg) must be flipped. Taking into account the irre- 
ducible corepresentation of the little group of /42m with 
the propagation vector k = (^,0, |), one can write the 
possible magnetoelectric coupling terms in the second de- 
gree of magnetic order parameters and the dielectric term 
in the Landau potential as: 

Fme = coEjPy+c'oElPy + c^EfPy+c^EfPy 

+ C Z 1 Z EfPy + C^EfPy + C\lElElPy 

\„xz rpx rpz p . xz rpx rpz ry . xz rpx rpz p 

+ c ll M/ + C 12 ^1 ^2 M/ + c 2l ^2 My 

+ c 22^2 E^Py + C 12-E 1 E^Py 

+ \{x- l P 2 x +Xy 1 P 2 v +X- z 1 P 2 z l (2) 

where the c coefficients denote phenomenological param- 
eters and Xxi Xy an d Xz are the components of the di- 
electric susceptibility. Since Mn magnetic moments have 
been experimentally observed in the xz plane, the terms 
containing Ei y and E2 y are omitted in the above equa- 
tion. The first and second terms denote exchange terms, 
the last term is the dielectric energy and the remaining 
terms are relativistic terms. The minimization of Eq.© 
gives a finite polarization along y: 

P y = - X y(coEi + c' E 2 2 + crEf + crEf 

1 ^zz rpz2 1 zz rpzZ , jxx rpx rpx 

H-c-l hj X + c 2 E 2 +c 12 E 1 E 2 

1 jxz rpx rpz , jxz rpx rpz , jxz rpx rpz 

+c 11 E 1 E 1 +c 12 E 1 E 2 +c 21 E 2 E 1 
+dSZE%E$ + cfiEtES), (3) 

whereas non-zero P x and P z are not obtained. This re- 
sult confirms Nenert's predictions: a spontaneous P y is 
evidently induced by the AFM coupling. However, at 
this point we cannot conclude whether the Heisenberg- 
exchange- or the relativistic-interaction terms is domi- 
nant in inducing P y . 
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FIG. 2: The upper panel shows Cu-3<i (a), Mn-3<i (b) and 
S-3p (c) partial DOSs of AFM Cu 2 MnSnS 4 and Cu-3d (d), 
Mn-3d (e) and Se-4p (f) partial DOSs of AFM Cu 2 MnSnSe 4 
(at U = eV). The lower panel shows the case of U = 4.0 
eV. The zero of the energy scale marks the Fermi energy. 



B. DFT results 

Figure [2] shows the partial density of states (DOS) pro- 
jected onto (a) Cu-3d, (b) Mn-3d and (c) S-3p states in 
AFM Cu 2 MnSnS 4 and onto (d) Cu-3d, (e) Mn-3d and (f) 
Se-4p states in AFM Cu2MnSnSe4. To discuss the elec- 
tronic structure, we focus on Cu2MnSnS4. The valence 
bands just below the Fermi level mainly come from Cu- 
3<i 10 states with small contribution from Mn-3d 5 states, 
while the bottom of conduction bands is mainly from 
Mn-3d states. We recall that, since Mn atoms are tetra- 
hedrally surrounded by S atoms, the five- fold degenerate 
3d states split into the three-fold £2 states and two-fold e 
states due to the crystal electric field (CEF). Addition- 
ally, the three-fold degenerated £2 states form bonding 
tb and anti-bonding t a states. The majority spin states 
are fully occupied whereas the minority spin states are 
fully unoccupied, corresponding to a high-spin d 5 config- 
uration (with 4.4 /iB magnetic moment inside the atomic 
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TABLE I: D FT- calculated P (/xC/cm 2 ) of AFM Cu 2 MnSnS 4 
and Cu2MnSnSe4 for U = and 4.0 eV, upon switching on 
and off SOC. When including SOC, all Mn spins are aligned 
along the x or z direction. 



sphere with 1.3 A radius). Cu-3<i states are weakly spin 
polarized with a small induced magnetic moment (equal 
to 0.015 /iB inside the atomic sphere with 1.3 A radius). 
We note, in fact, that the distance between Mn atomic 
sites is rather large, so that localized Mn-spins interact 
via Cu and S sites. Fig|2] (a) and (b) show that there is 
a small hybridization between Cu-3d and Mn-3d states, 
so that Cu-ions are slightly spin-polarized. As we'll show 
later, the small Cu spins caused by the Cu-Mn hybridiza- 
tion play an important role for ferroelectricity in these 
materials. 

Table U shows the DFT calculated P in the optimized 
structure, imposing the AFM configuration. Consistent 
with Landau theory results, only P y has a finite value. 
Both compounds show a rather small P compared to 
other materials with magnetically-induced ferroelectric 
polarization (comparingwith P equal to few /iC/cm 2 in 
E-type AFM HoMnOs[4], one or two orders of magni- 
tude smaller); this is likely due to the weaker "indirect" 
interaction (via Mn-S(Se)-Cu-S(Se)-Mn bond). 

The reason why Cu2MnSnSe4 shows a larger P than 
Cu2MnSnS4 has to be traced back to the fact that Mn- 
3d states can hybridize more with Se-4p states than with 
(more localized) S-3p states, in turn mediating a stronger 
interaction between Mn-3d and Cn-3d states. In the same 
respect, we note that the difference between U = and 
4.0 eV as far as P is concerned can be equivalent ly un- 
derstood in terms of Mn-3d and Cu-3d hybridization: the 
on-site Coulomb interaction pushes Mn-3d states deeper 
in energy (clearly shown in Fig|2]) and decreases the Mn- 
3d and Cu-3d hybridization. Therefore, P values at U = 
4.0 eV are smaller compared to the cases of U = eV. 

The important finding here is that the contribution 
to P due to spin-orbit coupling (SOC) is very small, as 
shown by the comparison of polarization values switching 
on or off the relativistic interactions (cf. Table [Ij). In this 
respect, we also note that the total energy decreases - by 
0.04 meV/Mn - when the spin direction is changed from 
the a to c axis, in agreement with experiments. 
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FIG. 3: (a) Total energy differences AE = E — Eo=qqo as a 
function of the relative angle at U — eV and U — 4.0 eV. 
(b) and (c) show P y calculated by the Berry phase method for 
the optimized atomic structure and for the atomic structure 
fixed to the case of = 90°, respectively. The lower panel 
shows the spin directions of Mn5, Mn6, Mn7 and Mn8 atoms 
(labels consistent with Fig[T]). 



We recall that our Landau theory analysis didn't clar- 
ify which magnetic term (i.e. Heisenberg exchange or rel- 
ativistic) mainly contributes to P; on the other hand, 
DFT results unambiguously show the relativistic term to 
give only a minor contribution. Therefore we can omit 
the terms depending on either E x or E z in Eq. ([3]), but 
keep the E 2 term so that the equation is simplified as: 



-XyicoEl + c'oEj) 

{-XyCo-XyC'ol^Si-Sj. 



(4) 



Here, P y depends only on the inner products of Mn 
spin vectors, therefore is described by the relative spin 
angle. Si- Si terms in the above equation are neglected 
because these terms disapper when calculating the polar- 
ization difference between FE and PE states. Progres- 
sive rotation of (£5, Sq, S?, S$) Mn spins can switch 
the ferroelectric polarization between +P and — P con- 
tinuously. Assuming the non-collinear (NC) spin con- 
figuration as Si = —S2 = —S3 = S4 = (0,0,1) and 
-S 5 = S 6 = -S 7 = S 8 = (sin (9,0, cos 6>), Eq.© be- 
comes 

P y = -x 2/ c (32cos6> + 32) - X?yCo(-32 cos (9 + 32). (5) 



Indeed, the P y oc cos 6 trend is obtained when we perform 
DFT calculations on the NC configuration by varying 
between and 180° with fully optimized atomic coor- 
dinates for each spin configuration. Figure [3] shows the 
calculated P^ p and the total energy difference AE{6) = 
E(0) — E{0 = 90°) along the adiabatic path as a function 
of the relative Mn-spin angle at Cu2MnSnSe4. We also 
show P^ p_Flx , calculated with atomic coordinates fixed 
to the case of 6 = 90°, so that there is no contribution to 
the electric polarization from atomic displacements. The 
similarity between polarization values for Figj3] (b) and 
(c) - including or not the ionic displacements - highlights 
stannite as a paradigmatic example of purely electronic 
ferroelectricity (see also discussion below). 

As expected from Landau theory and as pointed out 
above, the DFT-calculated Py P is fitted by a cosine 
curve, with coefficients ~Xy c o = Xy c o =0.0053 and 
0.0015 at U= and 4.0 eV, respectively. The func- 
tion AE shows convex and symmetrical behavior between 
positive and negative FE states. With increasing from 
= 0° (positive FE state), AE increases and reaches 
maximum value at = 90° (PE state) with the energy 
barrier AE(0)~ 0.68 and 0.25 meV/Mn at U= and 
4.0 eV, respectively. This energy barrier — coming from 
combination of JijSi-Sj term and P-related term — is 
much lower than in the HoMnOs case, where AE(0) = 
8 meV/Mn at U = eV.Q This is fully consistent with 
the weak Mn-Mn exchange-coupling constants involved. 



IV. MECHANISM OF FERROELECTRICITY 

According to the previous section, ferroelectric po- 
larization is induced by Heisenberg-exchange coupling 
between Mn spins. We also noted that the Heisen- 
berg exchange between collinear spins is more efficient 
in driving ferroelectric polarization than any spin-orbit 
Dzyaloshinskii-Moriya interaction [2lj: the order of mag- 
nitude of P here (tenth or hundredth of /iC/cm 2 ) is about 
one order of magnitude higher than that experimentally 
observed and calculated (hundredth or thousandths of 
fiC/cm 2 ) in spiral-based systems such as TbMnOa. The 
origin of ferroelectricity could be explained by analogy 
with the HoMnOs case, where O 2- ions move in or- 
der to reduce the magnetic coupling energy between Mn 
spins. Such "exchange-striction" mechanism has been 
referred to as "inverse Goodenough-Kanamori" (iGK) 
interaction. 0,[23 However, since here in Cu2MnSnS4 the 
Mn ions are far apart, the Mn spin interaction is too "in- 
direct" and it is difficult (and non- intuitive) to guess how 
S(Se) ions would move to reduce the magnetic coupling 
energy. Indeed, as shown in FigJS] (b) and (c), P^ p and 
P^ p_Flx have almost the same value: the ionic displace- 
ments are rather small and do not contribute much to 
P. 

Rather, to explain the "electronic" origin of the 
Heisenberg-driven polarization, we assume small-size 
spins (actually calculated as 0.015 //b) on Cu sites, which 
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(c) 

FIG. 4: Schematics of S(Se) sites tetrahedrally surrounded by 
one Mn, two Cu, and one Sn atoms in the 1 42m symmetry. 
Sn atoms do not have magnetic moment. The orange (large) 
and black (small) arrows indicate spins of Mn and Cu atoms, 
respectively. 



interact with Mn spins via "Cu-S(Se)-Mn" bonds. Con- 
sider S(Se) sites: from the chemical point of view, all 
of them are equivalent, since they are tetrahedrally sur- 
rounded by one Mn, two Cu and one non-magnetic Sn 
atoms in the symmetry 1 "42m. However, from the mag- 
netic point of view, due to the small magnetic moments 
at Cu sites, S(Se) ions can have two different magnetic 
environments at the FE state (6 = 0°): the magnetic 
moments of the surrounding Cu atoms can align either 
parallel or antiparallel, as shown in Fig0] (a) and (b). 
The anions are therefore split into two categories, they 
become inequivalent upon this peculiar AFM spin order- 
ing, breaking the symmetries and paving the way to fer- 
roelectric polarization. Since Mn spins interact via such 
magnetically inequivalent S(Se) ions, the charge densi- 
ties at Mn sites are modulated by Heisenberg-exchange 
to induce P. Similarly, charge densities at Cu sites with 
small magnetic moments caused by Mn-3<i and Cu-3<i hy- 
bridization can polarize. This is clearly observed in FigJSJ 
where we show the difference between FE and PE states: 
Mn and Cu charge densities are clearly non-symmetric 
along 6, i.e. P along y is induced. The results are con- 
sistent with the magnetically broken point group 2 at Mn 
and Cu sites (the point groups at Cu and Mn sites were 
"originally" - i.e. without magnetic ordering — 42m and 
—4, respectively). The "original" site symmetry at S(Se) 
sites is m, which is already a polar point group; therefore, 
a magnetically induced polar charge does not appear. We 
note that, in the PE state (0 = 90°), some Cu magnetic 
moments point along the x direction, corresponding to 
the NC Mn spin configuration. All S(Se) ions have then 
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FIG. 5: ab plane charge density difference between the FE and 
PE states in Cu2MnSnSe4: Ap(r) = />fe(V) — ppe(V). (a), (b) 
and (c) show the differences of charge density distributions at 
2=0, 0.125 and 0.25 planes (cf. FigUJa)), respectively. 



an equivalent magnetic environment, where half of the 
Cu magnetic moments align parallel to Mn and the other 
half perpendicular to it (Fig|H(c)); in this configuration, 
P is therefore not expected along any direction. 

V. SUMMARY 

We have investigated magnetically-induced ferroelec- 
tricity in AFM Cu2MnSnS4 and Cu2MnSnSe4 by means 



of Landau theory and DFT calculations. The combina- 
tion of both methods clearly highlights that the driving 
force of P is not the relativistic term, but the Heisen- 
berg exchange term. Microscopically, small localized Cu 
magnetic moments induced by Mn magnetic moments 
play an important role in inducing P. The Mn-Mn spin 
interaction via the magnetically inequivalent S(Se) ions 
produces a polar modulation of the charge density on 
Mn and Cu sites, finally inducing a net P with predom- 
inantly electronic contribution and negligible contribu- 
tion from ionic displacements. We therefore show that 
AFM Cu2MnSnS4 and Cu2MnSnSe4 are remarkable ex- 
amples of novel magnetically-induced ferroelectric mate- 
rials, with improper multiferroicity branching out into 
non-oxide and non-octahedral based systems. 
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